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Nonlinear State Observation Using 
/fo^-Filtering Riccati Design 

Konrad Rcif, Frank Sonnemann, and Rolf Lfnbehauen 



Abstract — la this paper, the authors propose an observer for 
continuous-time nonlinear systems. The observer gain is computed by 
a Riccati differential equation simitar to the extended Kalman filter. 
They prove that under certain conditions the proposed observer is an 
exponential observer by choosing an appropriate Lyapunov function. 
Furthermore, the authors explore some Important relations of the 
proposed observer to robust control theory and ifoc -filtering. To 
examine the practical usefulness of the proposed observer they applied 
it to an induction motor for the estimation of the rotor flux and the 
angular velocity. 

Index Terms — Extended Kalman filter, /foo- filtering, nonlinear ob- 
server. 



I. Introduction 

The basic model for the general nonlinear observation problem 
consists of a nonlinear plant and a nonlinear output map, at which 
neither system noise nor observation noise is considered. The problem 
is to estimate the states of the plant by observing Ihe output when 
the initial state of the plant is not known. In literature the design 
of nonlinear observers is a widespread area of current research (see 
[13] and [22] for a survey). Although many important results have 
been achieved, there are still possibilities for an improvement of the 
existing synthesis methods. 

A frequently applied method is the extended Kalman filter [5], used 
as a deterministic observer for nonlinear systems. The advantage is 
its easy handling with a simultaneous good performance; moreover, 
a full proof of its exponential stability has been given recently [16], 
[18]. However, a well-known difficulty arising through its application 
is the problem of divergence. The reason is that in computation of 
the observer gain only the linearized functions are considered, but 
not the whole non linearities. To obtain convergence, it is therefore, 
in general, necessary either to have a good initial guess so that the 
initial estimation error is sufficiently small or to permit functions 
only weakly nonlinear [16], [18]. 

Recently Hoc -control methods have attracted attention in robust 
control theory [4], [1 1]. It has been mentioned that the high robustness 
of die LQ-optimal control systems can be further improved by using 
i/ec -control [3]. In this paper we use a similar approach to improve 
the extended Kalman filter. Motivated by the results of robust control 
theory [9], [20] and tf^-ultering [2], [6], [10], [14]. The remaining 
nonlinear terms are treated as uncertainties and are implemented in 
the proposed observer design by a worst case valuation in such a 
way that these uncertainties can be tolerated. Consequently, better 
convergence properties can be expected; moreover, under certain 
conditions convergence can be achieved even independent of the 
initial estimation error. An alternative method for improving the 
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convergence properties of die extended Kalman filter has been 
proposed in [17]. 

The paper is organized as follows: in Section 11 we state the nec- 
essary preliminaries and introduce the proposed nonlinear observer. 
By choosing an appropriate Lyapunov function we show that the 
proposed observer is an exponential observer in the sense of [8]. 
Important relations to i/ec -filtering [2], [6], [10], [14] are explored in 
Section III. In Section IV we apply the proposed observer to estimate 
the rotor flux and the angular velocity of an induction motor (cf. [12], 
[15], [19]). Some conclusions are drawn in Section V. 

II. Proposed Observer 
Consider a nonlinear system represented by the equations 

= /(*(*)•*<*)) (1) 
y(t) = Cz(t) (2) 

where z(t) € R q is die state, x(t) € R p the input, y(t) € R m the 
output, and C a m x q matrix. The function /(•,■) and the input 
x{-) are assumed to be C 1 -functions. For the sake of simplicity we 
are restricted to linear time-invariant readout maps. The important 
case where some of the state variables can be measured directly is 
thereby treated. Furthermore, the linear time-invariant readout map is 
well suited to pointing out the main idea. Some remarks on the more 
general case with nonlinear readout maps are given in Section III. For 
the considerations following below we need a slightly more general 
definition of the exponential observer in [8], 

Definition: For a dynamic system given by (1) and (2) an expo- 
nential observer is a dynamic system 

i(0 = /(^W^'(0) + iSf(y(0^(0^(0^) (3) 

with an output injection g (%*.-,*) such that 

IK0- m\ < vlk(0)-i(0)||e-^ 

Vf. > 0, Vs(0) - i(0) e B fi (4) 

holds for some p,if,6 G R+ with B P = {t; € R <1 : ||-t> || < p}. 

Remark: The output injection #(■,■,■•') depends in general on the 
whole functions and x(-) and not only on the recent values z(t) 
and x(t). The reasons are given below. 

To calculate the output injection > •, *) we introduce a Riccati 
differential equation which is similar to the Riccati differential 
equation for the extended Kalman filter [5J. Using the abbreviation 

Ait) = ^(z{t).x{t)) (5) 

it reads 

P(t) = A{t)P(i) + P(i)A r (t) 

- P(t)(C T R- l C - X 2 I)P(t) + p 2 I (6) 

with some A, /* 6 Rr and a positive definite symmetric matrix R. 
We define die observer gain by 

K(t) = P{t-)C T R- 1 (7) 

and the output injection by 

y(v(t),z(-hx(-).t) = K(t)(v(t) - Cz(t)). (8) 

It is necessary to allow a general dependence on £(*) and x(-) and 
not only on the recent values z(t) and x(t) because the solution 
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of the Riccati differential equation (6) depends on past values of 

.4(0 = (df/dz)imMt)). 

An expansion of /(■,•) at z(t), .r(/) leads to 

f(z(t)Mt)) - f(m>*(t)) 

= A(t)(z(t) - i(0) + ip(s(t),z(t)**{t)) (9) 

where ip{z(t) 7 z (t).u-(t)) are the terms of second- and higher order 
in — zU)- Defining the estimation error by 

<(0 = *(0-i(*> (10) 
subtracting (3) from (1) and considering (8) and (9), we obtain 

at) = {A{t) - K(t)C)at) + v(z(t), f (*), Ht))- (id 

For an exponential observer this differential equation has an expo- 
nentially stable equilibrium point at zero; moreover, if (4) holds for 
all C(0) 6 Ii H the equilibrium point is globally exponentially stable 
(cf. pi, Sec. 5.1, p. 143]). 

To state the proof of the Main theorem given below we need the 
following preparation. 

Lemma: Consider real vectors z m z € R q and a* € R p , a symmetric 
q x q matrix IL and a function vK*, ■) with 



\\9{w)\\<k\\z-*\\ 



(12) 



for some n € R^. Choose \,fi £ R * such that k < A// holds. Then 
there is a 6 > 0 such that 

2(* - z) T U^(z, 5. ;r ) < (* - £) r [A 2 i + 0* 2 - S)UU] (z - z) 

(13) 

holds. 

Proof: From « < A/* it follows that there is a 0 < $ < A// 
such that 

n < A/* — 6 
holds. Inserting into (12) yields 

llv(*.M)ll < (A/i -MP -ill 

and therefore we have 

2(z - zfn?(z,z.x) < 2(A/i - s')\\U(z - z)\\\\z - (14) 
For considerations following below we evaluate the square 



(K) 



lin<*-*)||- A|| S 



= 2^1 + Q(*-i) r nn(*-i) 



||n(* - i)||||z - i|| + x\z - i) r (* - i) > o 



and use the fact that it takes only nonnegative values. Rearranging 
the terms yields 



2(V-A)||nu-i)||||U-i)|| 



nil + a 2 / 



(*-£). 

(15) 



Choosing 



6 = 



6__ 



from 6' < A// it follows that y > 6 /X and therefore 

Combining (14) and (1$) we obtain the desired result (13). □ 
With these preliminaries we are able to state Theorem 1 . 
Theomm J: Consider a nonlinear system with linear time-invariant 
readout map given by (1 ), (2) and an observer given by (3) and (6)-(8) 
with some A. // > 0 and a positive definite matrix R, such that the 
following conditions hold. 

1) For any solution £(-) of" the observer differential equation (3) 
the solution P(-) of the Riccati differential equation (6) is 
bounded via 



pi < R(*) < pi 



(16) 



for some p,p > 0. 
2) The nonlinear terms <p(z(t), z(t),x(t) ) in the differential equa- 
tion (1 1) for the estimation error are bounded via 



Mm*m,<t))\\<*\\*{t)-m\\ 



(17) 



with k < A/(. 

Then the considered observer is an exponential observer. Moreover, 
the differential equation (11) for the estimation error is globally 
exponentially stable at zero. 

Remarks: 

1) To ensure the global exponential stability of the differential 
equation (11) we require that (17) is satisfied in the whole 
state space. If this inequality holds only in a restricted area, only 
local exponential stability can be guaranteed. This means that in 
this case the initial estimation error may not be arbitrary large. 
Nevertheless, numerical simulations have shown that even if 
the stability is not globally satisfied, the proposed observer can 
tolerate much higher initial estimation errors than an observer 
based on the extended Kalman filter (see Section IV). 

2) Using similar techniques as in the proof following below, it 
can be shown that under certain conditions the continuous-time 
extended Kalman filter is an exponential observer [16], [18] too. 
The difference is that the equilibrium point of the estimation 
error is globally exponentially stable for the proposed observer 
but not for the extended Kalman filter. Tins means that the 
proposed observer can tolerate arbitrary large initial estimation 
errors, which is not the case for the extended Kalman filter. 

3) For the usual extended Kalman filter it can be shown that the 
solution of the Riccati differential equation is bounded via (16), 
if the matrices A(t) and C fulfill the uniform observability 
condition [7] or if the nonlinear system satisfies a certain 
delegability condition as stated in [1]. Unfortunately, these 
conditions cannot be applied to get similar results for the 
ifoc-fiitering-like Riccati differential equation (6). (The reason 
is that the term C T R~ i C - X' 2 I is in general not positive 
definite because C 1 R~*C does not have full rank if we treat 
nontrivial observation problems.) It is a well-known problem 
arising in -control as well as Hoc -filtering that the solutions 
of the corresponding algebraic Riccati equations and Riccati 
differential equations may lack being positive definite (cf. [3]). 

4) For the proof below as well as for the lemma stated above, 
the parameter ft does not play any important role. It has been 
introduced to have an additional degree of freedom during the 
implementation process. Especially for ensuring Condition 1, 
the additional design parameter /* has turned out to be very 
useful. Also see the following remark. 

5) The existence of a solution P(-) for the Riccati differential 
equation (6) which satisfies Condition 1 depends mainly on an 
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K{t) X {z{t)Mt)) 



Using 



$ — H / 



A(t) 



V (z{t).z(tU(t)) 



Fig. 1. Block diagram corresponding to the differential equation (27). 



w(t) 



[>(*)! 
[ml 



v(t) 



A{t)-K(t)C{t) 



linear part 



Fig. 2. Block diagram corresponding to the differential equation (31). The 
whole system consists of two interconnected subsystems. The dashed box 
contains the linear part of the system and the observer; Fit) and Hit) 
represent the higher order nonlinear terms. 



appropriate choice of A and //. As will be shown in Section 111, 
they are closely related to the attenuation constant 7 in H^- 
filtering problems. To find suitable values for A and // we 
employed a binary search algorithm, which is widely used to 
solve i/oo-control and /foe-filtering problems (see, e.g., [3]) 
and adapted it for our purpose (see Section IV). 
Proof: We consider the differential equation (11) for the es- 
timation error and prove its exponential stability by choosing the 
Lyapunov function 



v(CW»0 = C(*>n(0C(0 



(18) 



with 11(f) = P L (t). Because of Assumption 1 we have the 
following bounds for the Lyapunov function: 



i||C«)l| 2 < V{«t),t) < i||cC<)ll 2 



(19) 



which imply positive definiteness and decrescencc. The time deriva- 
tive of the Lyapunov function reads 

viath t) = < T mam) + c T (on(oc<o + c T ( wxw. 

Inserting <,'(?) according to the differential equation (11) yields 
V«(t),t) = C T (t)n(t)C(0 + [(A(t) - K(t)C)Cd) 

+ p(*(t).5(f),*(0)] r n(«K(0 
+C r it)n(t)i{A(i) - K(t)C)<;(i)+ v (z(t),z(t),x(t))]. 

Applying the lemma and considering (7) we get with a < X/i 

v(m,t)<<: r (tmt) + n(i)A(t) + A r (t)n(t) 

- 2C r R- 1 c + (p 2 - «)n(on(«) + a* /]<;(«)• 



ri(/) = -n(t)P(t)U(t) 



(20) 



(21) 



(22) 



and the Riccati differential equation (6) leads to 

v«{t),t) < < T (t)[-mw(t) - c T R~ 1 c\c(t). 

With the bounds (16) for P(t) we obtain the inequality 

v'(«o,o<-|riiaoii 2 

which implies that the time derivative of the Lyapunov function is 
negative definite. With (19) and (22) we satisfy the requirements to 
apply [21, Th. 62, Sec. 5.3, p. 173] and conclude that the differential 
equation (11) has a globally exponentially stable equilibrium point 
at zero. This includes that the coasidered observer is an exponential 
observer. □ 

in. Relations to #oo-Filtering and Robust Control Theory 

Considering the small gain theorem the proposed observer is 
closely related to #o©-filtering theory as will be shown below. Since 
these relations are more obvious for the case with nonlinear readout 
maps, we let a nonlinear continuous-time system be given by 

i(0 = /(;(0vK0) (23) 
y(t) = h(z{t)) (24) 

with a (^-function h(-). The linear readout map is thereby contained 
as a special case. For this nonlinear system we construct an observer 
according to (3) and (6)-(8)- Using the fact that the functions /(•,•) 
and h(-) are continuously diflerentiable, we expand them via (9) and 

h(z(t)) - h(z(t)) = C(t)(z(t) - z(t)) + *(*(*), i<0) (25) 

where C(t) is an m x q matrix given by 

C«)=fj<i(<)). (26) 

Defining the estimation error by (10) wc obtain from (23) and (3) 
using (8). (9), and (25) 

Cft) = (A(t) - K(t)C(t)X(t) + v»(*(0,i(0,*W) 

- K(i) x {z(t),m). (27) 

The nonlinear terms are assumed to satisfy (17) as well as 

\\ X (z(tym)\\<*\\m-m\\ (28) 

with some A > 0. Now we have a closer look at the differential 
equation (27) for the estimation error. Fig. 1 shows the corresponding 
block diagram. Defining 

1 



F(t) = 



i4K(t)\\ 2 



(29) 
(30) 



IICWII 2 

for C(0 i= 0 we have from (27) 

iit) = (A(t) - K(t)C(t) + DF(t) - K(t)H(t))C(t) (31) 
with D = \d as shown in Fig. 2. Moreover, we obtain 

[ml 

for every t > 0 which is an immediate consequence of (17), 
(28)-(30), and k < A/t. By this construction the nonlinearities are 
considered as time- varying uncertainties. For the stability of (31), we 
have the following theorem. 
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Fig. 3. Normalized angular velocity z-> (t) for the original system, the proposed observer, and the extended Kalman filter with a medium initial estimation error. 



Theorem 2: Consider (31) with D = /// and an observer gain 
A*(f) given by (6) and (7) with A./i > 0, and V? = i\ Let the 
following conditions hold. 

1) The Riccati differential equation (6) has a solution f J (*) which 
is bounded via 

V} < P{t) < pi 

for some p,p > 0. 

2) The matrices F{t) and H{t) are norm-bounded by 

-F(ty 



H(t) 



< A 



for every t > 0. 
Then the system (31) is exponentially stable. 

Proof: For the proof of this theorem we employ a variant of 
the small gain theorem. We recall the operator norms [10] for the 
two subsystems in Fig. 2 (e.g., [F J (t) H T (t)] r and the part inside 
the dashed box) 



||T wl .||oo = sup 

u-€L 2 

IMI 2 *o 



W|2 

1Mb 1 



||T vu .||oo = sup 
ll*'lla*0 



where || • || 2 denotes the usual La -norm for functions. According to 
[20, Corollary 1], the whole system in Fig. 2 is exponentially stable 
if the subsystem in the dashed box is exponentially stable and the 
small gain condition 

l|r VU! |U||r wt .-||oo < 1 (32) 
is satisfied. The operator norm HI^Hoo can be estimated easily to 

\m' 
[mi 



t>0 



< A. 



(33) 



Moreover, it can be seen that the subsystem inside the dashed box 
in Fig. 2 has the structure of a H. x -filter. Introducing 7 = 1/A and 



using R = J, D = ft I the Riccati differential equation (6) can be 
rewritten in the form 

P(t) = A{t)P(f) + P(t)A T {t) 

- P(t)C r (t)C{t)P(t) + P(f) f (<) + DD T (34) 

T 

and the observer gain (7) is given by 

K(t) = P(t)C T (t). (35) 

With (34) and (35) we are able to employ standard results in H^- 
filtering (cf. [2, Sec. 7.4, p. 301, Th. 7.6] or [6, Sec. 7.2.2, p. 267, 
Th. 7.2.1]) and conclude that the subsystem inside the dashed box 
in Fig. 2 is exponentially stable and the operator norm ||7Vw||oo is 
bounded by 

llilwlloc < 7. 

Together with (33) and 7 = 1/A we fulfill the small gain condition 
(32) and the proof is completed. □ 

IV. Application to an Example System 

To examine the practical usefulness of the proposed observer we 
apply it to a highly nonlinear example system, a symmetrical three- 
phase induction motor. The quantities to be estimated are the rotor 
flux and the angular velocity (see also [19] and the references cited 
therein). We introduce the state differential equations for an induction 
machine and do not mention too many details about the physical 
meanings of the considered quantities, which are contained in any 
standard book about electrical drives (see, e.g., [12, Sec. 10. 1] or 
[15, Sec. 3.3]). The state differential equations of a symmetrical 
three-phase induction machine read 

h(t) = kiZi{t) + zi{t)z 2 (t) + fe* 3 (0 + * 2 (t) (36) 
hit) = + £a(«) + faMO (37) 

i 3 (/) = k s z L (t) + k A zAt) + (a-i(t) - * 5 (0)*-i(t) (38) 
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Fig. 4. Lyapunov function V (£(£),<) for the proposed observer. 



i 4 (t) = k 3 z< 2 (t) - {xid) - Zo(t))zAt) + A^ 4 (0 P9) 
isft) = fes(si(0*-i(*) - Z2(t)z-4t)) + kexa{t). (40) 

All state variables are normalized. si(f), -22(f) and zs(t'), z.\(t) are 
the components of the stator and the rotor flux, respectively, in the 
plane perpendicular to the rotation axis. z^(t) is the angular velocity, 
xi(t) the frequency, and ^(O the amplitude of the stator voltage, 
respectively, and ara(t) is the torque of the load, ki , • * • , A*; are several 
parameters depending on the considered drive. 

The output are the normalized stator currents y L (t) and #2(0* 
leading to the output equations 



2/l(0 = *7SlW + fcsJ£ 3 (0 

2/2(0 = fe?s 2 (t)+ hz 4 (t) 



(41) 
(42) 



with parameters A7 and A«. The matrix ,4(0 can be calculated from 
(36)~(40), parameter £ omitted 
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and the output matrix is given by (41), (42) 



0 
A 7 



A* 
0 



0 

A* 



(43) 



(44) 



For the simulation we set ki = -0.186, A 2 = 0.178, A 3 = 0.225, 
A 4 = -0.234, h = -0.081, A 6 = -0.018, A 7 = 4.643, A* = 
-4.44$, a-i(0 = s{t). x 2 (t) = *(t), and .r 3 (0 = 0, where s(t) is 
the unit step function. These values are obtained, if in the original 
unnormalized system certain values for the physical quantities are 
chosen. To obtain appropriate values for A and /* we propose a binary 
search algorithm consisting of the following steps. 
1) Choose A,// > 0 and R - rl with some r > 0. 



2) Decrease r by taking the half- value until the solution P(-) 
of the Riccati differential equation (6) takes positive definite 
values. If no positive definite solution can be achieved, proceed 
with Step 3). 

3) Decrease yt by taking the half-value and initialize R = rl 
as in Step 1). Proceed with Step 2). If // gets smaller than a 
prescribed bound //u we let the algorithm abort. 

For the simulations we insert ,4(0 and C into the Riccati differen- 
tial equation (6) and choose X 2 - 0.7, /i 2 = 0.008, and R = 0.06/ 
according to the above stated search algorithm. The observer gain 
is given by (7) and die differential equation for the observer by (3) 
and (8). We solved the state differential equations numerically by the 
Runge-Kutta method choosing the initial conditions 

5(0) = [0.2 -0.6 -0.4 0.1 0.3] T (45) 

for the system to be observed and 

z(0) = [0.5 0.1 0.3 -0.2 4] T (46) 

for the observer. For the sake of comparison we also computed a usual 
extended Kalman filter [5] with the same weighting matrices ;* 2 J, R 
and the same initial condition i(0) by setting A = 0. /< 2 = 0.008, 
and R = 0.06/ in (6). 

Figs. 3 and 4 show die simulation results. As shown in Fig. 3 the 
estimation error approaches zero for the proposed observer but not 
for the extended Kalman filter. Further simulation results have shown 
that the proposed observer has often a larger domain of attraction than 
the extended Kalman filter, i.e., the initial estimation error may be 
larger. The Lyapnnov-function (18) is depicted in Fig. 4. Looking 
at the values for 0 < t < 1 we see that the Lyapunov function is 
increasing. This means that (22) is sufficient but not necessary for 
the error decay. 
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V. Conclusion 

The purpose of this paper is to propose an observer for nonlinear 
systems. This is done by modifying the well-known extended Kalman 
filter [5] for an effective treatment of the nonlinearities. We proved 
that under certain conditions this observer is an exponential observer 
[8] by choosing an appropriate Lyapunov function. For an alternative 
approach to the proposed observer we consider the nonlinearities 
as uncertainties and get a stabilization problem which is typical 
in robust control theory [9], [20]. This stabilization problem can 
be solved by ifo& -filtering techniques [2], [6], [10], [14] leading 
to the proposed modification of the extended Kalman filter. Using 
the proposed observer to estimate the rotor flux and the angular 
velocity of an induction motor we show that it can be easily applied 
to complex highly nonlinear systems. Numerical simulations show 
a good performance and an increased domain of convergence in 
comparison to the extended Kalman filter, i.e., we can tolerate 
larger initial estimation errors. Similar to ifoo -filtering problems, 
we have in general not a bounded positive definite solution of the 
corresponding Riccati differential equation, which docs exist for the 
extended Kalman filter if the system satisfies certain observability 
conditions [1]. Therefore, a careful choice of the constants in this 
Riccati differential equation is of particular interest. 
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Fractional-Order Systems and / J / A £)^-ControIlcrs 

Igor Podlubny 



Abstract — Dynamic systems of an arbitrary real order (fraction Hi- 
order systems) are considered. A concept of a fractional-order PI x D li - 
controller, involving fractional-order integrator and fractional-order dif- 
ferentiator, is proposed. The Laplace transform formula for a new 
function of the Mittag-Le filer- type made it possible to obtain explicit 
analytical expressions for the unit-step and unit-Impulse response of 
a linear fractional-order system with fractional-order controller both 
for the open and closed loop. An example demonstrating the use of 
the obtained formulas and the advantages of the proposed PI X D* 1 - 
controllers is given. 

Index Terms — Fractional-order controllers, fractional-order systems, 
fractional differential equations, Laplace transforms, transfer functions. 



I. Introduction 

Recently, several authors have considered mechanical systems de- 
scribed by fractional-order state equations [4], [5], [15], which means 
equations involving so-called fractional derivatives and integrals (for 
the introduction to this theory see [21]). 

There are also several recent applications in electricity. 

Le Mehaute and Crepy [14] have proposed a concept of a frac- 
tance — a new electrical circuit element, which has intermediate 
properties between resistance and capacitance. Such a device has been 
experimentally studied, for example, by Nakagawa and Sorimachi 
[18]. A circuit proposed by Oldham and Zoski [22] provides another 
example of a fractance. 

A new capacitor theory developed by Wcsterlund [28] is based on 
the use of fractional derivatives. 

New fractional derivative-based models are more adequate than the 
previously used integer-order models. This has been demonstrated, for 
instance, by Caputo [7], Nonnenmacher and Glockle [19], Friedrich 
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Development of a Servowriter for Magnetic Disk 

Storage Applications 

DANA H. BROWN and MUTHUTHAMBY SRI-JAYANTHA, member, ief.f. 



Abstract— The design, development, and application of a research 
servowriter* which is necessary to record magnetic position encoding 
schemes, are presented. The servowriter is designed to evaluate emerg- 
ing disk drive technologies. Submicron positioning accuracy of the disk 
actuator is achieved by integrating the state-of-the-art electronic, op- 
tic, electromechanical, and digital signal processing hardware. The 
servowriter also provides flexible and sophisticated signal encoding 
schemes through a programmable " pattern generator*' and has sev- 
eral modules of customized electronic hardware, each controlled by a 
supervisory computer. The position of the R/W (read/write) head and 
actuator assembly is measured by a laser interferometric system. The 
measured position is used to estimate the velocity of the R/W head, 
which is then used to implement a state- variable-based proportional- 
integral-derivative controller on the Intel (8086/8087) microprocessor 
chip set. Both systematic and random error sources associated with 
phase and amplitude encoding schemes are described. Environmental 
factors, such as seismic vibration and temperature change, which in- 
directly contribute to servowriting errors, are identified and steps to 
contain these error sources are described. 



1. Introduction 

THE HIGH -performance disk drives of today arc pre- 
cision rotating machines, with head positioning 
mechanisms capable of achieving submicron track fol- 
lowing accuracy. As the demand for track density in- 
creases, the ability to provide a higher resolution position 
signal will become mandatory. In order to maintain the 
magnetic R/W (read/write) head position along the center 
line of a desired track, a magnetic pattern that indicates 
the position of the R/W head is prerecorded on the disk 
surfaces during the manufacture of the disk drive. The 
device used to record the magnetic pattern is called a ser- 
vowriter fl, p. 719]. Since the magnetic pattern helps the 
servo control algorithms perform the track-following or 
track-accessing task, it is also referred to as a servo pat- 
tern. 

Of the several major steps involved in the manufacture 
of a disk file, the servowriting operation takes place just 
prior to sealing the head/disk assembly in a clean room 
environment. A servowriter in itself is a complex device 
that integrates state-of-the-art electronic, optic, and elec- 
tromechanical hardware with sophisticated signal encod- 
ing schemes to achieve the submicron accuracies that are 
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demanded in the industry. There are three broad catego- 
ries of position encoding techniques used by the disk drive 
industry: 

1) nonmagnetic position encoding (optical sensors 
mounted on the head actuator, mechanical detents 
for each track, stepper motors, etc); 

2) magnetic servo pattern based position encoding 
(amplitude encoded schemes [2] and phase encoded 
schemes [3]); 

3) combinations of the above (optical sensors with a 
single servo burst per track to remove low frequency 
track-following errors). 

The servo control of the IBM9332, for instance, is based 
on a phase encoded servo pattern [4], whereas the 
IBM33XX series is based on an amplitude encoded servo 
pattern. A short overview of servo patterns is given in the 
Appendix. 

This paper addresses the design, development, and ap- 
plication of a research servowriter that is targeted to 
achieve the magnetically based position encoding tech- 
nique in category (2). Following an overview of the ser- 
vowriter, its major elements are described briefly. The 
servowriter electronics are made up of several customized 
electronic hardware modules, including an Intel (8086/ 
8087) microprocessor chip set. Each module or subsys- 
tem is associated with a function that is required to sup- 
port the servowriting and testing process and is controlled 
by a supervisory computer (JBM-PC/AT). The pattern 
generator is software configurable and is controlled by the 
IBM-PC/AT. The IBM- PC /AT is menu driven and it is 
through this menu that the user controls or monitors the 
servowriting operation. More details are presented in |5|. 

The position of the R/W head is measured by a laser 
interferometric system [6]. The laser based position is 
used to estimate the state variables of the R/W head po- 
sitioning mechanism. These estimates are then used to 
implement a slate-variable-based proportional- integral - 
derivative (PID) controller, which positions the mecha- 
nism with zero steady -state error in the presence of bias 
force. Estimator/controller algorithms are implemented in 
8087-based floating-point arithmetic using assembly lan- 
guage. An example of a closed-loop transfer function of 
a servowriting mechanism is presented, and factors that 
would contribute to track-following error, which is det- 
rimental to the servowriter performance, are discussed. 
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Details related to clock pulse generation that is man- 
datory for locating the pattern generator pulses are de- 
scribed. The servo pattern is generated from 1024 x 32 
bit memory-based data, thus providing flexibility to study 
the performance of various encoding schemes. Finally, 
the error sources associated with encoding schemes, and 
environmental factors such as seismic vibration and hu- 
midity, which indirectly contribute to servowriting errors, 
are identified and steps to contain these error sources are 
described. 

IJ. Srrvowriter Description 
The servo patterns are written on one or more disk sur- 
faces, depending on the design of the file being written 
(i.e., dedicated versus sector servo). The servo patterns 
can be written with heads mounted on the disk drive ac- 
tuator, or a special set of heads can be mounted to the 
servowriter to perform this task. The servowriting process 
is a repetitive process of writing a track and moving the 
head to the next track. The process repeats itself until all 
tracks are written. The process of writing the servo pat- 
terns is relatively straightforward, but since it must be 
precisely done, the circuitry to do this can become quite 
complex. Fig. 1 shows the major elements of a servowri- 
ter, some of which are briefly described now. 

The heads must have a sensor mounted to them to in- 
dicate their radial position on the disk. Typically a laser 
interferometer is used to measure the head position. This 
sensor system is used in a closed-loop servo system to 
locate each track to be servowritten. Once the head is po- 
sitioned over a track, the pattern is written. However, the 
pattern must be synchronized to the disk rotational posi- 
tion so that adjacent concentric tracks line up according 
to a selected encoding rule. For this task, a sensor is used 
to supply fine resolution clock pulses that are synchro- 
nized to the angular position of the disk. This sensor is 
typically a stationary RMV head called a clock head. As 
the spindle motor speed varies due to random distur- 
bances, the clock pulses follow the speed changes accord- 
ingly. 

An encoder or pattern generator supplies the correct 
servo code to the heads for writing on the disk. This pat- 
tern generator follows the angular position clock as it 
writes each track so that the servo pattern on that track is 
aligned to the servo pattern that was written on the pre- 
vious adjacent track. The pattern generator is pro- 
grammed to change the code written on each track so that 
the servo signal on read back will indicate over which 
track the head is positioned. 

Each subsystem, together with the tasks performed by 
it, is coordinated by a central supervisory computer. A 
suitable architecture is used to share the computing, con- 
trol, and monitoring tasks required by a servowriter. 

III. Major Elements of a Servowriter 
A. Computer Architecture 

Fig. 2 shows the microprocessor hardware interconnec- 
tion necessary to achieve the servowriting operation. The 
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Fig. 2. Computer architecture. 

servowriter is controlled by an IBM PC/AT with external 
data output (DO) and data input (DI) registers. The con- 
trol program for this servowriter was written in IBM in- 
terpreter BASIC with some data manipulation programs 
(such as the pattern generator down-loader) written in 
8088 Assembly language. The memory in the PC/AT is 
640 kilobytes. However, the program to control the ser- 
vowriter is contained in a 64 kilobyte segment. Although 
interpreter BASIC is slow compared to Assembly lan- 
guage, 2000 servo half tracks can be written on a disk in 
well under 1 h. This is considered adequate for a research 
servowriter, but is not acceptable for a production ser- 
vowriter. The software would have to be optimized to re- 
duce that time by a factor of at least 5. 

B. Digital Position Control System 

The head is positioned by mounting a retro reflector to 
the actuator and tracking it using laser interferometry. A 
closed-loop digital servo feedback system is employed to 
hold the head position to within a few microinches of the 
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target position. A second-order actuator model with an 
augmented position integral state is used for controller de- 
sign. A moving coil actuator with a current-mode driver 
is thus modeled as follows: 
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where 

c s damping constant (N * s/m), 

k F actuator force constant ( N/ A ) , 

k s stiffness constant (N/m), 

k,. sensor gain constant ( laser counts /m ), 

m equivalent actuator mass (kgr), 

u control current commanded by servowriter pro- 
cessor (A), 

v laser measurement noise (equivalent laser count), 

w process noise vector (appropriate units), 

x T [position (m). integral (m • s), velocity ( m/s)], 

,x Rc( reference position (m), 

Z measured position (laser counts). 

The representative model defined by (I) is used to de- 
sign a PID-type controller using a pole-placement design 
technique. The introduction of process and measurement 
noise terms emphasizes that the system is not determin- 
istic, and that the tracking performance, therefore, can 
never be ideal (i.e., zero error). The unmodeled distur- 
bances act through the process noise term, h\ and mea- 
surement noise appears through the noise term, v. Seis- 
mic (or ground) vibration not only affects the accuracy of 
the model defined by (1). but also contributes to measure- 
ment noise by causing the laser optics to move relative to 
the actuator. Isolation of these sources contributes signif- 
icantly to the cost of a servowriter. 

The continuous model is discretized for controller/es- 
timator design. A sampling rate of 300 /*s is used to dis- 
cretize the model and is represented as follows ( ignoring 
noise terms, and setting x Rcr = 0): 

+ Tm, (3) 
z k = h T x k (4) 

where 
k sample index, 

h measurement matrix as defined in (2). 

Equations (3) and (4) are used to select the controller 
feedback gain constants. A third-order Butterworth pole 
pattern is used to locate the desired closed-loop poles of 
the position servo. 



Since the measured position (i.e., laser count) is not a 
continuous variable, a reduced order "current" estimator 
is designed to obtain the velocity estimate [7, p. 149]. In 
the design of the estimator, measured position is assumed 
to be exact, and the position (error) integral stale is ex- 
cluded from the estimator model. Therefore, the velocity 
estimator is of first order, and the pole-placement design 
yields a single gain corresponding to a desired estimator 
real pole location. The integral state is derived using the 
trapezoidal rule for integration. 

The controller/estimator gains and the system matrix 
elements are arranged together so that each laser position 
measurement is used to generate an appropriate electric 
current command in an explicit state variable form. The 
controller/estimator computations are realized in floating 
point arithmetic. Using a Hewlett Packard 3562A dy- 
namic signal analyzer, a closed-loop transfer function 
(CLTF) of an experimental system was obtained. Fig. 3 
shows the measured CLTF of an integrated disk/actuator 
system under servowriter control. The bandwidth of this 
system is limited by a 400-Hz resonance mode specific to 
the actuator used (in which the mounting of the retrore- 
flector accentuated the 400-Hz mode). The servowriter 
controller/estimator algorithm has adequate computa- 
tional speed to support a closed-loop system with much 
higher bandwidth if the actuator modes are not a problem 
from a stability standpoint. The presence of the integral 
state provides the required low-frequency disturbance re- 
jection. In order to suppress high-frequency disturbance, 
greater bandwidth with minimum CLTF peaking is rec- 
ommended. 

C System Timing and Clock Writer Design 

There arc different methods in practice today for writ- 
ing clock signals. The simplest of all methods is to let the 
clock writer run freely at a designed clock frequency and 
to dither the disk speed (because it can never be regulated 
perfectly at a constant speed) until closure (i.e., an inte- 
gral number of complete clock cycles around a track) is 
achieved. This task is automatically performed by the 
clock writer in 5 to 30 s. Hence, as defined in the next 
paragraph, the total number of clock cycles is first se- 
lected for a given problem, and through repeated write 
trials, the clock signal is recorded until satisfactory clo- 
sure is met. 

The total number of clock cycles per track is based on 
four parameters that have various system implications. 
They are defined as follows: 

A servo pattern sampling period (s), 

Wrpm spindle motor speed (rpm), 

/dock clock track frequency (MHz), 

Section integer number of sectors per revolution. 

The sampling time is specified to meet the performance 
requirements of the disk actuator servo. (Note that the ser- 
vowriter laser sampling rate is independent of the disk 
servo sampling rate.) The spindle speed is usually deter- 
mined from magnetics considerations and is set to a spec- 
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Fig. 3. Closed loop transfer function of an experimental system. 

itied value. The clock frequency should meet the R/W 
head frequency characteristics as well as the requirements 
imposed by the specific servo encoding scheme. The 
number of sectors is, however, constrained by the rela- 
tionship 

- 7T- (5) 

*-* ' * rpm 

The initial estimate of m xcxor will not necessarily be an 
integer. Therefore, a convenient integral number of sec- 
tors, m* ctor , corresponding to an exact, but slightly ad- 
justed, sampling rate A*, is selected as the first step. For 
instance, in order to drive the brushless dc spindle motor, 
an even number of sectors divisible by 12 (related to the 
number of motor poles) may be necessary. Based on the 
design requirements of the servo pattern, n xclot clock 
cycles at a frequency / clock arc to be generated between 
two consecutive sectors with a duration dictated by an ad- 
justed sampling rate A*. This requirement is met by the 
relationship 



= a* x y clock 



(6) 



In the case where « set(or of (6) does not result in an integer, 
it may be rounded to the nearest integral value, n* cct0I . 
Fractional adjustment of clock frequency (in megahertz) 
would typically satisfy this constraint. Hence, the total 
number of clock cycles per revolution of the disk spindle 
is 

" total ~ "sector ^ '"seaur (7) 

This redefines the spindle speed to be 

60 



^rptn ~ /clock * + 



X mZ 



(8) 



which, though not exact, is near the desired spindle speed. 
For example, a consistent set of parameters would be A* 
= 99.206 /zs, n^ t0T = 248 [-99.206 X 10' 6 x (2.5 X 
10 6 )], « tota , = 248 x 168 = 41 664 clock cycles per rev- 
olution, and N^ p - 3600.23 rpm, where the clock track 
writer is programmed to write n tota , (=41 664) clock 
cycles at a frequency of 2.5 MHz. This is done by turning 
on the write current to the clock head while writing 2.5- 
MHz data to it. In order to have exactly 41 664 clocks 
written under the clock head, a counter must register how 



many cycles have been written during each revolution of 
the disk. When exactly 41 664 clocks have been written 
for one revolution (achieved by dithering the speed of the 
motor slightly around Nrpm), the write current is shut off. 

Counting clock pulses during clock writing requires a 
reference mark to indicate precisely when the disk has 
turned through a single revolution. This is done by writ- 
ing an index mark with the servo head before clock writ- 
ing. The index decoder (see Fig. 1) monitors the read- 
back signal from the servo head and outputs a pulse each 
time the index mark is detected. This pulse (called index) 
initially starts the clock writer counter during clock writ- 
ing and resets it after each unsuccessful pass of clock 
writing. Clock track verification is done by simply read- 
ing back the written clock track and counting the number 
of clocks in one revolution. If disk defects damage one or 
more clock pulses, this will be reflected in the final count. 
An index mark is also encoded in the clock track to in- 
dicate the start of the clock track. Three cycles of 1.25- 
MHz frequency marker can be embedded in the clock 
track without seriously affecting the 2.5-MHz clock sig- 
nal. This index is used to synchronize the pattern gener- 
ators at the start of each servo track. 

The clock track resolution of 400 ns is too coarse for 
writing servo patterns, A type 2 phase-locked loop (PLL) 
that multiplies the clock track frequency by 16 and holds 
the output phase of this 40-MHz signal constant with re- 
spect to the clock track is employed. The PLL output ef- 
fectively subdivides the disk into 666 624 equally spaced 
time slots of 25 ns each. 

D. Motor Speed Control 

The speed control of the disk is critical during clock 
writing and clock read-back. If it varies too much during 
either clock writing or read-back, good closure will not 
be attained and the PLL that generates the system time 
base will go out of lock. This will make the writing of 
good servo patterns impossible. The speed of the disk is 
controlled by comparing motor Hall sensors (or, more ac- 
curately, the clock track signal) to a stable reference 
clock. For a research scrvowriter, the reference is derived 
from a high-precision (six to seven significant digits) fre- 
quency source. This is necessary for changing the number 
of clock cycles by a couple of counts per revolution. For 
5.25- or 3.5-in disk drives, the motors are small enough 
in power consumption that single IC motor controllers are 
readily available, and with the proper compensation, they 
can control the motor speed to about 0.01 percent of ref- 
erence speed. 

E. Pattern Generator 

The circuitry that writes the servo patterns on the tracks 
is called a pattern generator. It is usually programmable, 
but for servowriters with simple servo patterns it can be 
hardwired to write a limited number of servo codes. The 
pattern generator is a clock track counter that begins 
counting the 40-MHz PLL output clocks when the clock 
track index is received. At certain preprogrammed counts, 
the output of the pattern generator changes states to re- 
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verse the write current in the servo head. It is this se- 
quence of transitions that makes up the servo code. 

There are pattern generators of differing designs on the 
commercial test equipment market today. The pattern 
generator designed for a typical IBM servowriter is sim- 
ilar to some of the commercial ones except for instruction 
modifications to facilitate writing servo patterns. The out- 
put of the pattern generator is a shift register that clocks 
out the data with each 25-ns clock. The shift register can 
be loaded with 9-16 bits with each pattern generator in- 
struction and then shifted out at the system PLL clock rate 
(40 MHz). 

The pattern generator contains a block of memory, for 
example, 1024 X 32 bits. Each 32 bit word contains 3 
bits of instruction op code, 3 bits to indicate the shift reg- 
ister data length, 16 bits of data, and 10 bits of op code 
modifier (i.e., subroutine call address, branch addresses, 
or loop count information). Although the instructions are 
simple, they allow the writing of any servo pattern used 
to date and, in most cases, all tracks for an entire disk can 
easily be contained in the 1024 locations. 

F. Pattern Testing 

Writing servo patterns on a disk at exactly the correct 
locations does not guarantee a correctly written disk. It is 
possible that defects on the surface of the disk were pre- 
sent and that the servo information was not recorded at 
those points. The recording head may have had a shorted 
turn or other defect that prevented the servo pattern from 
being written on the disk also. 

The pattern test circuitry allows a read-back check to 
verify that a good servo pattern has been written. While 
the concepts and circuits for pattern testing are not diffi- 
cult, this part of the servowriter is troublesome. Errors 
are neither black nor white. They fall into a broad range 
of values, many landing on the threshold that the pattern 
test uses to evaluate the correctness of the pattern. The 
result is that a particular disk defect or servo error may 
pass on some tests and fail on others, or it may pass on a 
particular servowriter and fail on another. The pattern 
testing spans several subsystems and the details are pre- 
sented in [5]. 

IV. Error Sources 
The servowriter errors can be either random or system- 
atic and are the subject of the next section. 

Error Mechanisms in the Servo Pattern 

The phase encoded servo pattern looks at the average 
phase of many transitions with respect to a reference fre- 
quency (crystal). As a result, it is relatively immune to 
single transition defects and smalt amplitude variations. 
Ideally, the servo pattern phases should change continu- 
ously as the head position varies. Unfortunately the head 
gap geometry requires a discrete approximation. This ap- 
proximation to a sloping phase angle is a 5/2 ns discrete 
time shift for every t/2 track width of position change. 
Because of this discrete implementation, there can be rip- 



ple in the position signal generated by the demodulator as 
the head is displaced. Selecting a finer discrete resolution 
does not solve the problem of ripple completely. Mag- 
netic heads write tracks that are slightly wider than the 
width of the head gap. Each pass of the head while ser- 
vowriting leaves a noise fringe field between tracks. This 
fringe field is a fixed width that is a bigger percentage of 
each track width as track spacing dimensions shrink. If 
quarter tracks are written with 5/4 ns phase shifts, the 
linearity should get better (i.e., the ripple amplitude 
should decrease), but a larger percentage of the servo pat- 
tern would be made up of fringe field. Thus the improved 
linearity is more than negated by the noise induced. 

In amplitude servo schemes, the demodulator output 
versus the actual head position becomes nonlinear as the 
head moves away from the track center line. This nonltn- 
earity is very much dependent on the ratio of the head 
width to the track width. Thus the servo gain of the de- 
modulator output (millivolts/micrometers) varies signifi- 
cantly with head width. Some amplitude servo patterns 
use only a single transition in each position sample. If 
such a transition is damaged, though, the position signal 
for that servo sample would be incorrect. Longer burst 
servo patterns overcome this problem to a great extent by 
averaging several redundant transitions (as shown in Fig. 
6(b) in the Appendix). However, this redundancy cannot 
be carried too far because longer servo patterns begin to 
cause errors during high-speed seeks. They also consume 
the useful data space in files with sector servo architec- 
tures. 

1) Tangential Delay: Due to the skew of the head on 
the rotary actuator as shown in Fig. 4, the tops of the 
transitions on track N 4- 1 do not line up with the bottoms 
of the transitions on track N as shown in Fig. 4(a). How- 
ever, if the pattern written on track N is delayed slightly, 
they will line up as shown in Fig. 4(b). For such small 
adjustments, the 3-ns resolution is necessary. Besides 
head skew, nonradial motion of the actuator introduces a 
timing shift in the K, or tangential, direction, as shown in 
Fig. 5. The tangential delay circuit can be programmed 
to a new correction factor for each radial servo track to 
remove these systematic errors. There is still a random 
delay due to the uncontrolled vibration of the clock head 
(tocirtck) an ^ the servo head suspension (5y sen , 0 ). These 
two random components are time dependent, and are 
shown along the X and y direction for simplicity. The ef- 
fect of random vibration can be reduced by special treat- 
ment of the head suspensions or by selecting low jitter 
heads. Fig. 5 shows a generalized delay generating mech- 
anism relating the track spacing (fir), radius (r), disk ro- 
tational velocity (w), and slope of the actuator path 
(dy/dx). The relative time delay of a second track writ- 
ten next to an ideal first track (i.e., no head vibration) can 
be shown to be 

Relative Delay 

^(tyseivo - tow) + ^ SrJ 
: = [t X 5r] + (9) 
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Fig. 5. Effect of nonradial head motion on tangential delay. 

where r is the ideal delay per unit distance of r. There- 
fore, the first group of terms in the right-hand side of (9) 
represents the desired ideal time delay. However, in prac- 
tice the parameters defined in the second group of terms 
also contribute to the total delay. Observe that the delay 
parameters become more complicated when the actuator 
travels at higher velocity, as in the case of a seek. 

2) Spindle and Head Suspension Jitter: As pointed out 
above, track to track timing alignment is essential. Me- 
chanical movement of the spindle center line as it rotates 
will cause the disks to shift the track center line contin- 
ually under the servowriter heads, even if the access sys- 
tem holds the heads perfectly stationary. Radial nonre- 
peatable spindle motion can cause the servo pattern being 
written to squeeze toward one of the adjacent tracks. Tan- 
gential nonrepeatable spindle motion can cause the data 
to be written out of alignment with the adjacent track, 
even if the timing signals were to occur at the proper 
times. Head suspension jitter can cause symptoms similar 
to spindle jitter. The head suspension jitter comes from 
both the clock head and the servowriting head. The jitter 
primarily occurs at right angles to the head suspension's 



main axis, as the stiffness of the suspension is weakest in 
that direction. In-line heads reduce tangetial jitter, 

3) Vibration Isolation and Mechanical Alignment: In 
order to minimize file vibration during servowriting, the 
file is isolated from external vibrations by mounting it on 
a large granite base that floats on air isolation mounts. 
Ideally, the file should be clamped at three points, so that 
the base casting is not distorted while servowriting. In 
practice, bending modes of the casting sometimes can oc- 
cur with just three clamps holding the file down and these 
vibrations can be written into the servo pattern in the form 
of track positioning errors. Thus some experimentation 
will usually have to be done to determine the best clamp- 
ing arrangement for a particular disk drive. 

The alignment of the disk and head stack for sector ser- 
vowritten files is critical. During file operation, all heads 
on an actuator should be centered over the same track 
number on each surface. In addition, the servo patterns 
on all surfaces should pass under the heads at the same 
time. At the time of servowriting, this alignment is veri- 
fied. However, if the file is clamped incorrectly in the 
servowriter or in the final product configuration, the ac- 
tuator and disk spindle will tilt, resulting in radial mis- 
alignment and timing shifts from head to head. Sector 
servo files are designed to allow for some radial and tim- 
ing misalignment, but typically these numbers are on the 
order of 5 /*m and 500 ns, respectively. 

4) Environment: The key to successful servowriting is 
cleanliness and the precise placement of servo informa- 
tion. For this reason, the environment of the servowriter 
is extremely important. A servowriter should be in a clean 
room so that exposed files will not become contaminated 
with dirt particles. Even in a clean room, though, dirt will 
accumulate on any surface (although minimal). A simple 
cover should be placed over each open file in the clean 
room while it is watting for an assembly operation or ser- 
vowriting. During the servowriting operation the cover 
can be removed temporarily. Besides keeping dirt off of 
the disks and out of the head-disk assembly (HDA) en- 
closure before it is sealed, magnetic particles should be 
purged from the clean room. Temperature is critical in 
precision placement of the tracks. The file dimensions will 
shift as the temperature changes, so it is best to keep the 
temperature between 68 and 72 °F to minimize these me- 
chanical changes and to facilitate operator comfort. Hu- 
midity should be between 40 and 60 percent. Above 60 
percent, the humidity causes the servowriter to rust. Be- 
low 40 percent, static electricity begins to cause prob- 
lems. 

V. Conclusions 
Servowriters may seem complex on the surface because 
many interleaved subsystems are involved. However, each 
subsystem is easy to understand. Because of the building 
block approach, the servowriter is very logically con- 
structed. It is a precision machine that should be treated 
with care. Time and money spent on the development of 
a servowriter usually are recovered in a less expensive 
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product. However, a servowriter cannot be designed and 
developed independently of the hie. The servowriter's 
mechanical clamping and access systems are directly tied 
to the mechanical design of the file. Likewise, a file may 
behave very poorly after servowriting if its design has not 
been thoroughly verified on the servowriter. However, no 
matter how much planning and testing are performed, in- 
evitable adjustments to the file and servowriter do occur 
during development. They usually come in the areas of 
vibration isolation, file alignment, and actuator position- 
ing. For this reason, mechanical design considerations in 
the file and servowriter should be given high priority. En- 
hancing either the software or the electronics alone will 
not compensate for mechanical design deficiencies. The 
environment causes problems that are not directly trace- 
able. Magnetic particle contamination may appear once 
the product has left the servowriter. Errors in servowrit- 
ing induced by temperature shifts or head crashes due to 
accumulation of debris never directly point to the cause, 
thus requiring anxious search for solutions, usually later 
during high-volume production. With care, it can be an 
extremely powerful analysis tool, and without it, a never- 
ending source of problems. 

Appendix: Servo Patterns 
The phase encoded servo pattern (or phase modulated 
servo) works on the principle of varying the relative phase 
of two bursts of squares waves (both at the same fre- 
quency) for each successive track. As a head crosses 
tracks, the phase shift between burst 1 and burst 2 for each 
track changes. This phase shift becomes an indicator of 
track position. In Fig. 6(a) the vertical dashed line rep- 
resents the boundary between burst 1 and burst 2. When 
the head is positioned over track N + I, there will be no 
phase shift between burst I and burst 2. When the head 
moves to track N + 2, there is a 90° phase shift between 
burst 1 and burst 2. The phase shift varies quite linearly 
with position for all head locations between tracks N + I 
and N + 2. Note that the pattern repeats every four tracks 
so that the same phase shift results when the head is over 
track A/, N + 4, N + 8, • • • . The benefits of a phase 
encoded pattern are that the servo signal frequency does 
not vary, the pattern is linear over four tracks, and single 
bit disk defects usually do not have a significant effect on 
the servo position signal, since multiple cycles are aver- 
aged. 

The amplitude servo scheme is more commonly em- 
ployed today. It has good linearity over track centers, but 
becomes nonlinear as the head moves off track center by 
more than 35 percent of the track width. This scheme is 
usually easier to demodulate and to servowrite. Ampli- 
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tude servo relies on the principle that signal amplitude 
read by a head varies linearly as the head moves off track 
center. Referring to Fig. 6(b), the read-back signal am- 
plitude from the head will be at its maximum during burst 
A when the head is positioned over track Vv\ N + 2, TV + 
4, • • • . As the head moves off track, the amplitude of 
the burst A signal decreases in proportion to the amount 
of off-track position, and burst B amplitude begins to in- 
crease from 0. At tracks AT + \.N + 3 t N + 5, • • , the 
signal from burst B is at its maximum and burst A is zero. 
Like the phase encoded servo, this scheme will tolerate 
small disk defects since many cycles arc used in the burst. 
However, amplitude servo is not as linear as a phase en- 
coded servo and there are two widely separated frequency 
components in the pattern— dc and the burst frequency. 
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Abstract: 

In industrial motor drive systems, a shaft torsional vibration is often generated when a 
motor and a load are connected with a flexible shaft. This paper treats the vibration 
suppression control of such a system. In this paper, a control system design method 
using linear matrix inequality (LMI), which is a tool for control design that replaces or 
complements Lyapunov-Riccati equations, is provided and the speed controller for an 

induction motor by LMI is proposed. In the H w speed controller, weights are used to 
satisfy tracking and disturbance rejection. Experimental results show the validity of the 
proposed speed controller by LMI, and this controller is compared with the state 
feedback controller 
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Abstract - In the industrial motor drive system, a 
shaft torsional vibration is often generated when a 
motor and a load are connected with a flexible shaft. 
This paper treats the vibration suppression control of 
such a system. In thi fa papfei ' , a contr ol systeiri design 
method using Linear Matrix Inequality (LMI), which 
is a tool for control design that replaces or 
complements Lyapunov-Riccati equati ons, is provided 
and the H» speed controller lor an induction motor 
by LMI is proposed. In the H« speed controller, 
weights are used to satisfy tracking and disturbance 
rejection. Experimental results show the validity of 
the proposed H« speed controller by LMI, and this 
controller is compared with the state feedback 
controller. 

1. Introduction 

In the industrial motor drive systems such as 
industrial robots, rolling mills, elevators and so on, 
the shaft torsional vibration is often generated when 
a motor and a load are connected with a flexible 
shaft. Such a system is called as the two-mass 
resonant system and this paper treats the vibration 
suppression control of this system. 

To suppress the shaft torsional vibration, the 
resonant ratio control, the state feedback control, the 
LQG control and the H« control are studiedti-ej It 
is often impossible to detect the load speed, the 
shaft torsional torque and the load disturbance 
torque in the real systems. Recently, it is studied 
the new method which estimates unknown state 
variables by using the reduced order observer and 
feed back these state variables by using the pole 
placement design method^] But this design method 
doesn't satisfy the fast speed response property and 
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the attenuation property of load disturbance at the 
same time. 

In recent studies, there is the H«> control which 
satisfies the fast speed response property and the 
attenuation property of load disturbance at the same 
time by using the weights. The H~ control is often 
applied in the industrial systems. 

In recent, because it is possible to represent the 
control specifications and the robust stable conditons 
by Linear Matrix Inequality (LMI), the control 
problem can be transformed into the convex 
problem.^] LMIs have gained much attention as a 
tool for control design that replaces or complements 
Lyapunov-Riccati equations. 

In this paper, the controller design method by 
LMI is provided and the H» speed controller of an 
induction motor by LMI is proposed. The validity of 
the proposed controller is comfirmed by 
experimental results and this controller is compared 
with the conventional state feedback controller using 
the pole placement method.^) 

2. The two-mass resonant system 

The two-mass resonant system which has a 
motor and a load connected with a flexible shaft is 
shown in fig. 1. The state equation of the 
two-mass resonant system is as follows 
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Fig.l The plant model of the two-mass 
resonant system 



Jm 

0 
0 



Jl 
0 



(1) 



where, 

Jm. Jl : the motor inertia, the load inertia 

B M t &l ' the motor viscous damping coefficient, 

the load viscous damping coefficient 
K S Ff : the shaft stifmess 
Wm > <*>i : the motor speed, the load speed 
T M , T S h f Tl : the motor torque, the shaft 

torsional torque, the load disturbance torque 

The goals for the control of the two-mass 
resonant system are stated in the following forms, 

(1) Suppress the shaft torsional vibration. 

(2) Reject the effect of the load disturbance torque. 

(3) For the change of the speed reference, the load 

speed follows it quickly without overshoot. 



23 



PLANT P(») 
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Fig. 2 The generalized plant of the H*> control 

T L is the disturbance torque and n is the speed 
command. W d is weight for T L and W n is weight 
for n. If weights are defined as follows, 



L <- Wn U Wn\ 



Wu 



(3) 



then the state space equation for the generalized 
plant and the ratio of output to input, G , are as 
follows. 



(4) 
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(5) 



3. The H~ generalized plant 

Fig. 2 is the H«> generalized plant to apply the H 
» control theory to the two-mass system. In fig. 2, 
iv y iv t are inputs for disturbance and command, 
respectively, and z u z u are outputs. z x is as 
follows. 
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(2) 



Above G is obtained by choice of W d , W n , W u 
and the H~ controller K is solved by Matlab. 

4. The control system design by LMI 

In this chapter, we represent G of the 
generalized plant by LMI. Convexity of LMI and 
LMI representation of Riccati inequality for the 
control design are given and convex problems 
involving LMI are treated. 
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Riccati inequality applying to control system 
design is as follows. 

A T S+ SA + SBR~ l B T S + Q<Q , /?>0 (6) 
In Eq. (6), A and B are 



A = 



A, 0 BjCw 
0 A m 0 
0 0 i4 w 



0 BJDiw Bp 
B m 0 0 
0 B w 0 



(7) 



in real two-mass system and Q, R, S can be 
sellected by free choice. Eq. (6) can be transformed 
into 



<* r 2 7-F r KF>0; V: — W T W> 0 (10) 
where W : scaling matrix, r : level 

following problems can be obtained. 

1) CFP : when matrix F and level r>0 are 
given, problem solving scaling matrix W( or 

v ); 

3 V: = W T W>0 such that y 2 K- F r 7F>0 

2) CMP • when matrix F and scaling matrix W 
are given, problem minimizing level y\ 

minimize r> 0 such that y 2 V - F T VF >0 

We can get W( or V ) in case of minimum y 
by iteration. 



-AS- SA - Q SB] _ 



B'S 



> 0 (8) 



by using Schur corriplementuo] And Eq. (8) can be 
converted to following general form of LMI. 



Fie) = F 0 + 2 ei F t 



(9) 



where 



e, : vector 



-A T S f -SiA SiB 
B T Sj 0 



Convex problems for solving F(e) > 0 are 
remarkably two problems as follows. 

1) Convex Feasibility Problem (CFP) 
find ee?? m such that He) > 0 

2) Convex Minimization Problem (CMP) 
minimize c T e subject to F\e) > 0 , eeK m 

where ceH m is given. 

And because of scaled matrix norm condition, 
II ^FPF" 1 !! < r, can be transformed into eq. (10), 

II wfvt 1 \\ < r » vr~ T f 7 * w^wf frr -1 < r 2 / 



5. EL» speed control of the two-mass system 
by LMI 



We design the H~ controller of the two-mass 
system by LMI. In fig. 3, the weighting function 
W u is selected as W u = 10 ~ 6 which does not 
influence the design of controller. The choice of 
weighting functions is very critical because the 
control performances are determined by weighting 
functions. The weighting function W d is for the 
closed-loop transfer function from the disturbance 
torque 7/, to Z\. In order to achieve the design 
specifications such as vibration suppession and 
disturbance rejection, | W d | " J is chosen as 

closed-loop gain plot of the one-mass system. The 
weighting function W n is for the closed-loop 

transfer function from the disturbance torque n to 
Z\ , where n corresponds to the speed command or 

noise. The closed- loop transfer function from n to 
wu is a complementary sensitivity function, so 
W n is chosen as a high pass filter considering the 

robust stability. Then Wd and W n are chosen as 

follows. 
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Fig. 3 Bode diagram of weights 
( ver : [dB], hor : [ 0) ] ) 

d 7 d (s+KT 3 )(s+10 5 ) ' 

W = y (J + 5Q) (11) 

Wn Yn (5+10 5 ) U 

The H« controller K(s) can be obtained by using 
Matlab after aajusting the parameters in weighting 
functions as 7 d = 5, w d = 40 , r „ = 500. 
Weighting functions are shown in fig. 3. The H« 
controller K(s) is given by 

K(s) -4.921 xlO'x 

^+lxifl 5 )(^n.6±?8Q)(s+5.41) 

(5 + 5.6X 10 5 )(s + 3.59* 10 4 )(s+9±;33.52)(s+0.0011) 

(12) 

6. Experimental system 

The configuration of the experimental system for 
the two-mass resonant system is shown in fig. 4 
There is a flexible shaft ( 1 : 800 [mm], <t> : 13 
[mm] ) between a motor and a load generator. This 
system is used variable inertia to change the inertia 
of a motor or a load generator. For this paper, it is 
only used variable inertia for a load generator. The 
mechanical constants of the two-mass resonant 
system are listed on Table 1. 




Fig. 4 The experimental system for the two-mass 
resonant system 



Table 1 The mechanical constants of the two-mass 
resonant system 



1m 


Jl 




0.008 Dcgml 


0.08 [kgm*] 


50.527 [Nm/rad] 



7. Experimental results 

Experimental results show the validity of the 
proposed controller and this controller is compared 
with the state feedback controller. The inertia ratio 
is 0.008 : 0.08 ( ] M : J L ) in this paper. 
Experiments are composed of speed response 
experiment and load disturbance response 
experiment. 

7,1. Speed response experiment 

The reference speed of the load is 200[rpm] as a 
step input at 0 [sec]. Fig. 5, 6 are experimental 
results of step response under the state feedback 
control. Fig. 5 shows load speed and motor speed 
and fig. 6 shows shaft torsional torque. The speed 
transient time from 0 to 200 [rpm] is about 0.7 
[sec]. If the transient time becomes shorter, 
overshoot of the load speed is larger. In the 
transient time, maximum torque is 6 [Nmj. 
Fig. 7 ( 8 are experimental results for control by 
LM. Because the speed transient time from 0 to 
200 [rpm] is about 05 [sec], the transient time of 
proposed LMI controller has much less than that of 
state feedback controller. This is because of 
adjusting the weighting function for speed response. 
In the transient time, maximum torque is 9 [Nm]. 
But this maximum torque does not influence the 
two-mass system 
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Fig. 5 State feedback controller 
150 [rpml/div, 0.5 sec/div 
(upper) Real motor speed, (down) Real load speed 



Fig. 9 State feedback controller 
150 [rpml/div, 0.5 sec/div 
(upper) Real motor speed, (down) Real load speed 




Fig. 6 State feedback controller 
3 [Nml/div, 0.5 sec/div 
Real shaft torsional torque 




it s itimV u t todt a i / 



Fig. 10 State feedback controller 
3 [Nml/div, 0.5 sec/div 
Real shaft torsional torque 
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Fig. 7 LMI controller 
150 [rpml/div, 0.5 sec/div 
(upper) Real motor speed, (down) Real load speed 



Fig. 11 LMI controller 
150 [rpm]/div, 0.5 sec/div 
(upper) Real motor speed, (down) Real load speed 




Fig. 8 LMI controller 
3 [Nml/div, 0.5 sec/div 
Real shaft torsional torque 



Fig. 12 LMI controller 
3 [Nml/div, 0.5 sec/div 
Real shaft torsional torque 
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7.2. Load disturbance response experiment 

A resistor load which can be connected with DC 
generator is used as a load. By using a resistor 
load the load disturbancs torque is increased from 0 
to 4.5 [Nm] at 200 [rpm]. 4.5 [Nm] is 37% of rated 
torque. Fig. 9, 10 are experimental results for the 
state feedback control. Load speed dropped 17% of 
200 [rpm] and recovered after 1,2 [sec]. This is 
because the inertia ratio is big. If the inertia ratio 
increases, the recovery time is much longer. 

Fig. 11, 12 are experimental results for control by 
LML The load speed dropped 5% of 200 [rpm] and 
recovered after 0.6 [sec]. Then, the speed drop of 
the proposed LMI controller has much less than that 
of state feedback controller. The recovery time of 
the proposed LMI controller has much less than that 
of the state feedback controller. 



8. Conclusion 



This paper treated the speed control of the 
induction motor system which has a motor and a 
load connected with a flexible shaft. The H«> speed 
control of an induction motor with two-mass 
resonant system by LMI is proposed. In the 
proposed control, the H« control, which satisfies the 
fast speed response property and the attenuation 
property of load disturbance at the same time by 
using the weights, is used and is repesented by 
LMI which is a tool for control design that replaces 
or complements Lyapunov-Riccati equations. 

The experiments show the excellent performance 
of the proposed H« speed control by LMI compared 
with the state feedback control. 
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